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ABSTRACT 

We have used widely spaced in time Hubble Space Telescope images to determine tangential ve- 
locities of features associated with outflows from young stars. These observations were supplemented 
by groundbased telescope spectroscopy and from the resultant radial velocities, space velocities were 
determined for many outflows. Numerous new moving features were found and grouped into known 
and newly assigned Herbig Haro objects. It was found that stellar outflow is highly discontinuous, as 
frequently is the case, with long term gaps of a few hundred years and that these outflow periods are 
marked by staccato bursts over periods of about ten years. Although this has been observed in other 
regions, the Orion Nebula Cluster presents the richest display of this property. Most of the large scale 
Herbig Haro objects in the brightest part of the Orion Nebula appear to originate from a small region 
northeast of the strong Orion- S radio and infrared sources. With the possible exception of HH 203, 
we are not able to identify specific stellar sources, but do identify candidate sources for several other 
bright Herbig Haro objects. We find that there are optical features in the BN-KL region that can be 
related to the known large scale outflow that originates there. We find additional evidence for this 
outflow originating 500-1000 years ago. 

Subject headings: Herbig Haro Objects — Star Formation — Galactic Nebulae:individual(Orion Nebula, 
NGC1976) 



1. INTRODUCTION 

The long operational life of the Hubble Space Tele- 
scope (HST) has presented the opportunity to obtain 
and compare emission-line images of the central (Huy- 
gens) region of the Orion Nebula over a sufficiently long 
time-base to determine highly accurate motions in the 
plane of the sky. These motions enable us to study the 
outflow that accompanies early phases of star formation 
and the interaction of this out-flowing plasma with the 
ambient nebular material. An introduction to the Orion 
Nebula and its associa ted cluster is presented in several 
recent review articles flO'DellpOOH |Muench et aT1|2008 " 
O'Dell et al.|[2008 ). 



In this paper we present the results of the determi- 
nation of new tangential and radial velocities and their 
interpretation. We will show that a number of new stel- 
lar outflows (Herbig Haro objects, HH) have been found, 
that some previously designated HH objects are not what 
was originally argued, and that these motions provide us 
with important new information about star formation in 
the visually obscured BN-KL and Orion-S centers. 

2. OBSERVATIONS 

Our observational material consists of both new spec- 
troscopic and imaging data. The spectroscopic observa- 
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1 Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc., under NASA Contract No. NAS 5-26555. 

2 Based on observations at the San Pedro Martir Observatory 
operated by the Universidad Nacional Autonoma de Mexico. 



tions were made at Observatorio Astronomico Nacional 
en San Pedro Martir operated by the Universidad Na- 
cional Autono ma de Mexico using t he MEZCAL echelle 
spectrograph ( Meaburn et al.||2003|) . The new imaging 
observations were made with the HST's WFPC2 and 
ACS cameras as part of program GO 10921. The data 
were processed with the IRAF package of software. 3 



2.1. Spectroscopy 

The Huygens region of the Orion Nebula has been the 
subject of many spectroscopic studies. Most of these 
did not have velocity resolutions to determine with suf- 
ficient accuracy the radial velocities of features moving 
about the speed of sound (about 10 km s _1 ) and the 
studies with this necessary spectral resolution usually 
treated only certain sections of the nebula. More com- 
plete high resolution mapp ing of the Huygens region have 
been made recently (|Doi et al.|2004l|Garcia-Diaz fc Hen- 
|ney|2007 Garcia-Diaz et al. 2008) in the strongest lines. 
Those multi-slit-setting maps were made with primarily 
north-south oriented slits, so that the spatial resolution 
was excellent along the slit length, but poorer in right 
ascension because the spacing was usually wider than 
the characteristic "seeing" . In this study, we draw upon 
some of those earlier spectra, but our primary target was 
an extremely detailed east-west region centered on the 
Orion-S area and new observations were necessary. 



3 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. under cooperative agreement with the 
National Science foundation. 
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Fig. 1. — TheselST^x 26. 4" samples of a WFPC2 image (O'Dell & Wong 1996) have superimposed the locations of the 15 slit spectra 
obtained. The declination shifts in arcseconds relative to 159-35U are shown tor each setting. Below each image is a contour plot of the 
strong (6583 A) [N II] line beginning at 2000 counts and increasing in steps of 3000 counts. By setting the initial contour at 2000 counts one 
sees the velocity changes in the nebula more clearly, at the expense of the visibility of the faint high velocity features. The short horizontal 
line in each contour image is the position for a heliocentric velocity of 20 km s — . Full resolution versions of all of the figures in this paper 
are included with the electronic edition of this journal article. 



These new observations were made over the period 11- 
14 January 2007. The instrument configuration and per- 
formance was the same as in o ur earlier study of the 
Ring Nebula ( jO'Dell et al.|2007| ). The spectrograph was 
employed in a single spectral order mode by isolating 
the order containing the Ha 6563 A and [N II] 6548 A 
and 6583 A lines with a narrow-band interference filter. 
The 70 /im wide entrance slit projected 0.9" onto the 
plane of the sky. We processed the well exposed sample 
of 187". The 24 /am pixels of the Site 1024x1024 CCD 
detector were double-sampled along the slit (0.62" per 



pixel). The spectrum was read out in single pixels along 
the dispersion direction, giving Full Width at Half Max- 
imum (FWHM) of the Thorium- Argon comparison lamp 
emission lines of about 7.0 km s~ . The atmospheric see- 
ing varied from one to one and one half arcseconds during 
the exposures. A series of east- west oriented slit settings 
were made using offsets from the bright star 159-350= J W 
499 (in the p osition based syst em introduced in [O 'Dell 



fc Wen| (|1994D or the catalog of | Jones fc Walker | JTS55JJ 



and at intervals of about 1" at locations shown in Figure 
[I] We will use the position based system of designation 
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throughout this paper, adding decimal places when an 
accuracy of greater than about 1" is needed. Exposures 
of 600 s were made with characteristic peak signals of 
30,000 counts in Ha. Transparency variations were cor- 
rected by normalizing each slit's signal against a similar 
region of a flux calibrated ([O'Dell fc Doi| |1999 ) image of 
this object flO'Dell fc Won g 1996) made with the HST's 
WFPC2, but gaussian blurred to resemble ground-based 
image quality. Figure [I] shows contour plots of the re- 
sulting spectra. 

Because of the high velocity resolution and excellent 
signal level, it was possible to search for faint high ve- 
locity features. The results of this search are shown in 
Figure |2j which superimposes the velocities of high ve- 
locity components upon the slit settings and a WFPC2 
image. For completen ess, this data also shows the results 
from an earlier study (Doi et aT||2002). 



early and later images. This i s the same technique used 



in our earlier investigations (Bally et al. 2000; Doi et 



al. 



2002 O'Dell & Doi 2003) using data witha shorter 



2.2. Imaging 

Our imaging program was planned to allow a homoge- 
nous set of data to be obtained, which was intended to 
allow the best possible comparison of earlier and later 
images. In planning the observations it was recognized 
that the high spatial resolution of the ACS-WFC (0.05"/ 
pixel in the drizzled images) was better than that of the 
WFPC2 (0.09967 pixel), but the latter instrument al- 
lows a longer time base since it was installed earlier. The 
planning became more complex as the ACS-WFC failed 
and WFPC2 observations were substituted where possi- 
ble. 

2.2.1. New Observations 

In this paper we report on the analysis of new images 
in four areas of the Orion Nebula, supplemented by simi- 
lar investigation of pairs of earlier observations. The four 
areas were selected to coincide as much as possible with 
earlier images of what were thought to be the most in- 
teresting and relevant objects. The characteristics of the 
new observations are summarized in Table [TJ The total 
exposure times presented were obtained in triple expo- 
sures, which allowed removal of cosmic ray induced spots 
on the images. On-the-fly routinely processed data were 
the starting point for our study. 

2.2.2. Earlier Observations 

We drew on a large body of earlier observations made 
with the WFPC2 and the ACS-WFC. Their character- 
istics are presented in Table [2] This listing includes not 
only the images matched with the new observations, but 
also images that were reprocessed and analyzed as part of 
this program. In addition to these long time-base ACS- 
WFC images, we also employed archived F658N ACS- 
WFC images obtained through program GO 10246 with 
Massimo Robberto as principal investigator. 

2.2.3. Processing the Image Data 

The first step in comparison of the first and second 
epoch images was to align the two. In each case we took 
the first epoch images as the standard and shifted the 
second epoch images. Since the individual CCD detec- 
tors of both cameras can change with respect to their 
neighbors, alignment and comparison was only made us- 
ing reference stars falling into a single CCD in both the 



time interval. The accuracy of the alignment was typ- 
ically about 0.2 pixels, which corresponds to about the 
3 km s~ internal di spersion of velocities of the Orion 
Nebu la Cluster stars ( Jones fc Walker|1985 Furesz et al. 



2008). Because the first and second epoch fields of view 



were not exactly on the same part of the sky, the region 
useful for comparisons is more irregular and smaller than 
the second epoch images. 

After alignment the signal of the two images was nor- 
malized and a search was made for moving objects by 
taking a ratio of the first and second epoch images. In 
the case of identical images, the ratio image would be 
a smooth value of unity and in the case where objects 
are moving there would be a pattern with a dark leading 
edge and a bright trailing edge. In most cases motion of 
individua l features was m easured using the least squares 
method (Hartiga n et al.| |2001) employed previously, or 
when the objects were large or the regions complex, sim- 
ply measuring their positions directly. The measurement 
uncertainties ranged from 0.1 to 0.4 pixels, depending 
on the nature of the object being measured. The un- 
certainty in the direction of the motion is less well de- 
fined, but can be 30° for very small motions and only 
a few degrees for the large motions. The results of the 
measurement of discrete objects are presented in Table [4] 
(Appendix B) and Table [5] (Appendix C). However, there 
were many large-scale moving features not measured in 
either fashion, but the reality of the motion in the ratio 
images is unquestionable. 

Figure 3 and Figure 4 show ratio images for our Visit 
1 and Visit 2 images for the WFPC2 F656N and F658N 
filters, supplemented by images near HH 203 and HH 
204 from a comparison of archived GO 5193 and GO 
8121 images which have been smoothed by median filter- 
ing. The difference in time between the first and second 
epoch images was 12.76 years for Visit 1, 12.49 years for 
the Visit 2 F631 and F673N images, 13.71 years for the 
Visit 2 F502N, F547M, F656N, and F658N images, and 
6.25 years for the HH 203/204 images. At our adopted 
distance of 440 pc (c.f. Appendix A), these intervals cor- 
respond to WFPC2 velocity scales of 16.3, 16.6, 15.2, and 
33.2 km s _1 per WFPC2 pixel respectively. Ratio images 
in F502N were created and used for finding rapidly mov- 
ing objects but are not shown here as all moving features 
in this filter are also seen in the F656N image, whereas 
many of the F656N moving features are not seen in the 
F502N image. Figure 5 shows the same tangential veloc- 
ity image superimposed on a WFPC2 color image. 

The detailed coding of the superimposed data for Fig- 
ure 3, Figure 4, and Figure 5 are described here, rather 
than in the figure captions. The markings on the outer 
borders show the Declination south of —5° and the Right 
Ascension east of 5 h 35 m . Filled circles indicate the po- 
sitions of stars mentioned in the discussion and those 
within the dashed outline are found in the near infrared 



catalog of Hillenbrand & Carpenter (2000). Arrows indi- 
cate the direction and magnitude of the tangential veloc- 
ities, with the arrow length corresponding to the motion 
in 300 years. The colors of the arrows are added only to 
indicate local groupings, the same color being used mul- 
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Fig. 2. — This is the same WFPC2 image as Figurel but showing all of the slit settings. The labels on the short axis are the Declination 
values south of 5° while the labels on the long axis are the Right Ascension values east of 5 h 32 m (both epoch 2000). The -I, -II, and -III 
refer to features within the multiple eastward moving shocks that form HH 529. The letters J show the locations of spectral features that 
must be microjets since they show a range of velocities starting at the system velocity, blue indicates Ha and red [N II]. The numbers 
indicate the heliocentric radial velocities and all of them are negative, that is, the minus signs have are not shown. The green numbers 
indicate high velocity features in Ha and/or [OIII] found in|Doi et al.|(|2004|). The other numbers show results from our slit spectra. 



TABLE 1 

Program GO 10921 Data Used in this Study 



Visit 


Target 


Camera 


Date 


Modified Julian Date 


Filter 


Exposure (seconds) 


1 


Orion-S 


WFPC2 


2007-12-25 


54459.3 


F502N 


300 


- 


- 


- 


- 


- 


F547M 


90 


- 


- 


- 


- 


- 


F631N 


600 


- 


- 


- 


- 


- 


F656N 


180 


- 


- 


- 


- 


- 


F658N 


300 


- 


- 


- 


- 


- 


F673N 


480 


2 


HH 528 


WFPC2 


2007-09-16 


54359.8 


F502N 


400 


- 


- 


- 


- 


- 


F547M 


120 


- 


- 


- 


- 


- 


F631N 


600 


- 


- 


- 


- 


- 


F656N 


240 


- 


- 


- 


- 


- 


F658N 


400 


- 


- 


- 


- 


- 


F673N 


540 


5 


HH 502 


ACS-WFC 


2006-11-22 


54061.4 


F502N 


1300 


- 


- 


- 


- 


- 


F658N 


945 


- 


- 


- 


- 


- 


F660N 


1000 


6 


HH 540 


WFPC2 


2007-11-05 


54409.7 


F502N 


1000 


- 


- 


- 


- 


- 


F547M 


80 


- 


- 


- 


- 


- 


F656N 


600 


- 


- 


- 


- 


- 


F658N 


800 


- 


- 


- 


- 


- 


F673N 


800 



TABLE 2 

Archived Observations Used in this Study 



Program 



Target 



Camera Date Modified Julian Date Filter Exposure (seconds) 



GO 5193 Orion-Southeast WFPC2 1993-12-29 



GO 5469 HST1,HST10,LV3 WFPC2 1995-03-21 



GO 8121 



GO 8121 



GO 9460 



HH 203-204 WFPC2 2000-03-20 



BN-KL 



WFPC2 2000-09-13 



HH 502 ACS-WFC 2002-08-07 



49350.i 



49797.5 



51633.1 



51800.9 



52493.7 



F502N 


800 


F547M 


200 


F656N 


700 


F658N 


1200 


F502N 


300 


F547M 


90 


F631N 


300 


F656N 


180 


F658N 


300 


F673N 


300 


F502N 


780 


F656N 


780 


F658N 


1800 


F673N 


1800 


F502N 


780 


F656N 


780 


F658N 


1800 


F673N 


1800 


F502N 


2000 


F658N 


1000 


F660N 


1000 
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Fig. 3. — The determination of tangential motions for the Huygens region are shown in this 213 // x219 // mosaic of F656N ratio images 
(first/second), with north up. See the text for details of the labeling. The nearly square central dashed region is shown enlarged in Figure 
7 and the elongated dashed region is the spectroscopy field in Figure 1 and Figure 2. 
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Fig. 4. — The same as for Figure 3, except this shows the F658N ratio images. 



tiple times for different groups. The dashed elongated 
rectangular outline indicates the region of radial velocity 
results shown in Figu re 2. Open red c ircles indicate the 
position of H2 knots (Stan keet al.|2QQ 2 J) and red lin es H2 
features seen in Subaru images pvait'u et al.||2QQQ ). The 



long irregular blue lines delineate collimatea high veloc- 



ity flow as indicated by He 1 1 0830 A emission (Takami et 
aLl |2002| [Henney et al.||2007| ). The st rong dark red/ blue 
contoured lines indicate CO outflow ( |Zapata et al.|2005| ) 
and the pastel colore d pink/light-blue co ntoured lines in- 
dicate SiO outflow ([Zapata et al.|[2006|). Op en sq uares 



in dicate the positions of H 2 maser sources (Gaume et 
al. 1998). Point sources within the dashed outline are 
coded by the shortest wavelength of their detection, with 
filled white squar es indicating the positions of radio only 



rilled wnite squares indicating tne positions 01 radio only 
visible sources ( |Zapata et al.|2004a|b 2005 ) , red squares 
the positions 01 sources seen onfy~m th e mid infrared 



(Smith et al. 2004 Robberto et al. 2005), and filled or- 



a nge circles positions of stars in the near infrared catalog 



of [Hillenbrand fc Carpenter (2000). A few motions out- 
side o f our field of view ha ve been added with the c oding 
OD= (|0 , Dell fc Doi||2003j), DOH =j |Dpi et al.|[2002| , and 
G-DH=(Carcia-Diaz & Henney 2007]). The region most 



likely to contain the sources of the high velocity optical 
outflows is shown as a dark dashed ellipse. 

Figure 6 shows the Visit 5 region of overlap with the 
GO 9460 field near HH 502. In this case the time in- 
terval was 4.29 years and the velocity scale 24.3 km s _1 
per ACS-WFC pixel. It has had the faint fine features 
enhanced by division by a third order polyn omial fit to 
the ba ckground. The nomenclature is that of |Bally et al. 
( [2006 ). 

3. DISCUSSION 

This rich body of new data allows us to draw conclu- 
sions about many aspects of gas-flow and star formation 
in the Orion Nebula Cluster (ONC) and its secondary 
star formation regions in the BN-KL and Orion-S re- 
gions. In this discussion we will first consider outflows 
associated with the Orion-S region (both as individual 
systems and as a possible global unit). We then treat 
what the new data tell us about the BN-KL region. The 
multiple outflows to the south of the Huygens region and 
centered on HH 502 are discussed third. The discus- 
sion concludes with consideration of individual proplyds 
and outflow systems. Throughout this paper we will as- 
sume that the velocity of the host Orion Molecular Cloud 
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Fig. 5. — The same as for Figure 3, except the background image is a section of the WFPC2 mosaic image (|Q'Dell fc Wong 1996). 



(OMC) is ^=25.8 km s" 1 JQ'Dell et al.|2Q08| > and that 
of the ONC is 25.6 km s" 1 flFuresz et al.|2008| ). Velocities 
relative to the OMC will be designated as VIomc • We will 
continue to assume that the distance is 440 pc, so that 
one-tenth parsec subtends 46.9". Position Angles (PA) 
will be expressed in the ordinary way, counter-clockwise 
from north and #, the angle of flow with respect to the 
plane of the sky, is positive in the direction of the ob- 
server. We draw heavily on publis hed slit-spectroscopy 
maps of high er ionization (|Doi et al . 2004) and lower ion- 
ization lines Garcia-Diaz fc Henney| ( |2007p (henceforth 
DOH04 and GD07), and where we do not have new tan- 
genti al velocities , on|D oi et al.| (|2002[), henceforth DOH, 

^01 



or 



O'Dell fc Doll ( [2005] ), hencefortn~6D03. 



3.1. Outflow from the Orion-S Region Near-Surface 
Sources 

There is by now a rich literature about large-scale 
(greater than 0.1 pc) outflows probably originating from 
the Orion-S region, the idea going back to the symmet- 
ric forms of HH 202 and HH 203+HH 204 being pointed 
in opposite directions across Orion-S, then developing 
as more outflows (HH 269, HH 529, HH 528, HH 625) 



were found. Henney et al. (2007), henceforth HO07, have 
most recently summarized the history and characteristics 
of these outflows, including their possible connection to 
much larger apparent shocks outside of the Huygens re- 
gion, so that it is not necessary to discuss each object 
anew, beyond summarizing how the new observations 
have altered our understanding. Our new observations 
do give valuable new information about the exact region 
of the sources of these outflows. 

Once again we find the well known result that all of the 
large-scale outflows are blue-shifted with respect to the 
OMC. When their current motions are projected back- 
wards, these projected lines cross in the Orion-S regions, 
suggesting that they have a common origin in one or 
more objects located there. The results for these objects 
are summarized in Table [3] The region that probably 
contains the sources of the high velocity optical outflows 
is shown as an ellipse in Figure 3, Figure 4, and Figure 5. 
This region is a more refined estimate of the location and 
proper ties of the OOS (Op tical Outflow Source) identi- 
fied in |OT)ellfe Doi| ( |2003| ) . 



HH 202. — HH 202 is at the apex of a large parabolic 
high ionization feature, as shown in Figure 5. The variety 
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Fig. 6. This 187"xl78" inverted brightness ACS-WFC F658N image of the region around HH 502 has PA=-23.4° as up. Motions 
over a period of 100 years are shown. The dashed boxes show the regions depicted in larger scale in Figure 10 and Figure 11. The vertical 
strip with many small dark features is a zone of poor removal of cosmic ray events. 



TABLE 3 

Motions of Large-Scale HH Objects Originating near Orion-S 



Designation 


V T 


source 


Vq 


source 


Vector Velocity 


e pa 


(degrees) 


Timescale (years) 


HH 202-South 


59 


GO 10921 


-38 


DOH04 


89 


48 


326 


2,600 


HH 202-North 


45 


GO 10921 


-21 


DOH04 


67 


47 


336 


3,600 


HH 202-NW 


31 


GO 10921 


-18 


DOH04 


55 


56 


310 


5,600 


HH 202-WNW 


38 


GO 10921 


-20 


DOH04 


61 


52 


298 


5,100 


HH 203 


77 


DOH02 


-46 


DOH04 


77 


43 


139 


3,400 


HH 204 


92 


DOH02 


-18 


DOH04 


92 


26 


137 


3,400 


HH 269 


50 


GO 10921 


-46 


GO 10921 


88 


55 


275 


1,800 


HH 269-Ram 


18 


GO 10921 


-25 


GO 10921 


54 


71 


290 


1,500 


HH 528-Base 


22 


GO 10921 


17 


HO07 


24 


22 


178 


5,700 


HH 528-Cap 


24 


GO 10921 


17 


HO07 


24 


22 


159 


5,500 


HH 529-III 


40 


GO 10921 


-30 


GO 10921 


69 


54 


119 


2,000 


HH 529-11 


92 


GO 10921 


-31 


GO 10921 


108 


32 


91 


780 


HH 529-1 


54 


GO 10921 


-30 


GO 10921 


78 


46 


90 


980 


HH 529-Ram 


59 


GO 10921 


-44 


GO10921 


84 


45 


95 


410 



Note. — All velocities are in km s and Timescales are calculated from the tangential motions as if the point of 
origin was at 144-352. GO 10921 indicates data from this study. 
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of ionization stages indicates that although it forms as a 
collimated high velocity outflow it also contains neutral 
material. This can be due to the flow impacting neu- 
tral material or that the flow has compressed the ambi- 
ent ionized gas to the degree that it traps the ionization 
front. If the latter, this would require that the mate- 
rial lies within the main cavity of the nebula rather than 
in the foreground veil. The white arrows in this region 
in Figure 5 indicate tangential motions in the HH 202- 
South brightest part of the object and the red arrows 
the HH 202-North part. The yellow arrow indicates flow 
in a shock feature ahead of HH 202-North (designated as 
HH 202-NW) and the blue and orange arrows two shocks 
even further ahead (designated as HH 202- WNW). The 
collimated flow probably driving these shocks has been 
measured by DOH04 to have radial velocities of about 
Vomc — — 40 km s _1 to —60 km s _1 and the d etails of 
this flow ar e best shown in the He I 10830 line ( |Takami| 
et al. 2002), which is seen in good contrast against the 



background emission from the nebula because the line is 
usually quite optically thick except where the emitting 
materially is Doppler shifted off the core of the line. 

HH 203. — HH 203 is a well defined bow shock with 
low ionization features at its tip. It is apparently driven 
by a high velocity, high ionization jet that emerges into 
the zone ionized by ^ 1 Ori C or the nearer (in the plane 
of the sky) # 2 Ori A (which dominates the ionization of 
some of the proplyds in this area) about half the distance 
from the point of origin. The jet radial velocity is Vq = 
-28 krns" 1 . 

HH 204.— HH 204 is nearly at the same PA as HH 
203, but shows marked differences. Its apex is a floccu- 
lent structure that includes low ionization features even 
though the body of the enclosed parabolic form has ex- 
tended [O III] emission, indicating that it is forming in 
ionized gas. No driving jet is seen. It should be noted 
that HO07 found a broader shock front extending beyond 
even HH 204 and at a PA midway between HH 204 and 
the terminus of HH 528, as if there had been three events 
of collimated flow in about this direction. 



HH 269. — HH 269 is composed of a series of well defined, 
wide- ionization ra nge, shocks oriented almost due west 
( Walter et al.|1995| ). In this study, only features near 116- 
345 were measured (the so-called HH 269-East shock) but 
both in tangential and radial velocity. A 5" narrow but 
wiggly low ionization feature was also measured at its 
western terminus and is designated here as HH 269- Ram. 



HH 528. — HH 528 shows a broad band of irregular, low 
ionization, structures headed southeast from the Orion- 
S region. It is composed of two regions, the base or 
"jet" (as designated in HO07) and the southeast-most 
bow shock (again as designated in HO07). The term 
"jet" only loosely applies as this feature is much broader 
than any collimated flow producing the other HH ob- 
jects in this study and one sees individual bow shocks 
within the object. Falling b ack on the original discovery 
work notation for HH 528 ( |Bally et al"]|2000D where it 
was noted that the entire object is shaped like a mush- 
room, we will here designate the "jet" as the "base" and 
the end feature as the "cap". HH 528 may in fact rep- 
resent two outflows, the "base" object defined by broad 



irregular structure, oriented towards PA=155° and con- 
taining features moving towards PA=178° and a second 
outflow represented by the "cap" , oriented very approxi- 
mately towards 147° and containing features moving to- 
wards 159°. In both cases the objects may be passing 
through the low-ionization portion of the nebula, rather 
than beomg within the main body of photoionized gas. 
A mediating argument against division of HH 528 is that 
the velocities and spatial orientation of the two compo- 
nents are very similar. 

There is a peculiar complex of features around 
5:35:17.6 —5:24:54 (all coordinates in this paper will be 
in epoch 2000). Parts of this complex seem related to 
the HH 528 flow, such as the 177-454 shock, which is 
measured in [S II] to be moving at 37 km s _1 towards 
PA=199°. However, there are other features that seem 
to be moving at velocities of order 130 km s _1 towards 
PA = 100°. Most puzzling is a feature just to the south- 
east of 177-454, where a naive interpretation of the [N II] 
ratio map would imply an oppositely directed motion of 
130 km s~ towards PA = 280°. However, careful in- 
spection of the individual images in multiple lines sug- 
gests that this feature is also moving towards PA= 100°. 
The reason for the discrepant behavior in [N II] is that in 
this filter the feature seems to be a moving dark feature 
rather than a moving bright feature, which reverses the 
sense of motion in the ratio image. In addition, the mor- 
phology of the feature seems to be evolving, becoming 
darker in the second epoch. It is much larger and more 
diffuse than the shock features associated with HH 528 
and at this point remains unexplained. 

HH 529. — HH 529 is composed of a series of shocks 
oriented towards the east and moving in that direction. 
Here we break it down into designated shocks HH 529-III, 
HH 529-11, and HH 529-1 proceeding from east to west (as 
shown in Figure 2-Figure 5). There is a series of small, 
rapidly moving features which may represent the driving, 
nearly collimated outflow that we designate as HH 529- 
Ram. There may be an HH 529-0 shock which is not 
obvious on images, but reveals itself through high radial 
velocity features a few arcseconds north and south of the 
HH 529-Ram group (cf. Figure 2). The HH 529-Ram 
features begin at a Right Ascension of 5:35:14.91 and 
end at 5:35:15.6 while the putative HH 529-0 shock is at 
5:35:15.5, i.e. near the west end of the Ram features, and 
its velocity Vr - 
V T = -44 km s 
Bow shocks HH 529-III, HH 529-11, and HH 529-1 are 
high ionization only, indicating that they fall within the 
zone of ionization by ^ 1 Ori C, a conclusion supported by 
the detail ed spectroscopic an alysis of HH 529-III+HH 



—48 km s is comparable to the Rams' 
, which also argues for an association. 



529-11 by flBlagrave etlI1[2006| . However, the HH 529- 
Ram feature is seen in both high and low ionization lines, 
indicating that this outflow is moving from a low into a 
high ionization region. 

3.2. Outflow from the Orion-S Deeply Imbedded Sources 
3.2.1. Molecular Outflows 

The brightest infrared and radio sources lie to the 
southwest of the putative source or sources producing the 
multiple large-scale HH objects associated with Orion-S. 
Recent radio observations reveal two well defined molec- 
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ular outflo ws in this southwest part of Orion-S. CO ob- 
servations (Zapata et al. 2005) detail a bipolar outflow 
with the blueshifted component towards PA=305° and 
apparently arising from the infrared star 136-40 0, which 
has associated H2O masers ( |Gaume et al.||1998J ). A sec- 
ond study in SiO reveals a bipolar molecular outflow 
that apparently arises from the infrared star 136 -355, 



with its blues hifted component towards PA=284° (Zap- 
ata et al. 12006). These features are shown in Figures 3-5. 



There are several spots of SiO emission in the strong ra- 
dio sources 137-408 (CS 3) and 134-411 (FIR 4), but no 
obvious directed outflow. 136-400 and 136-355 have no 
optical counterparts and they must lie behind a signifi- 
cant amount of obscuring material, so one would expect 
to find optical components to the molecular outflows as- 
sociated with only the blueshifted components. 

The SiO blueshifted component coming from 136-355 
intersects with the H H 269 series of shock s, which have an 
orientation of 273° ( Walter et al.||1995 ) , an orien tation 
strengthened by the presence of long H2 filaments ( |Kaifu| 
|et al.|[2Q00| ). The feature for which we have measured 
the tangential velocity lies upon a projection of the SiO 
blueshifted outflow, but the direction of motion of these 
features (PA=275°) more nearly agrees with the overall 
orientation of HH 269 and a line projected back to the 
HH 269-Ram feature. 

3.2.2. The Optical HH 625 Object Associated with the High 
Velocity CO Outflow 

There is a well defined series of optical and infrared 
features associated with the CO blueshifted component. 
One sees H2 features extendin g towa rds 136-400 and 
there are a series of H2 knots (Stank e et al.||2002[ ) ex- 
tending out to 204" (0.44 pc). There is a peculiar, low 
ionization, optical feature that has been designated as 
HH 625 and we measured four features in it with an av- 
erage motion of Vt=34 km s _1 towards PA=296°. We 



have examined Ha and [N II] spectra of this region ( Doi 



et al. 2004) and determined that V© = —8 km s~ . This 
means that this optical feature (which must be formed 
in a region of lower obscuration than applies closer to 
136-400) has a velocity vector of 48 km s~ inclined 41° 
towards the observer with respect to the plane of the sky. 
The CO observations trace the blueshifted component to 
Vq = —36 km s~ (^omc = — 62 km s~ ), which means 
that the optical feature probably represents mass-loaded 
material moving more slowly than the original outflow. 
However, the small size of the H2 knots indicates that 
the flow that drives them has remained very narrow and 
collimated. 

3.3. Outflows from the Central Region 

In Figure 7 we show the central portion of the Huy- 
gens region that includes Orion-S. This incorporates all 
of the area where the axes of symmetry of the flows (or- 
dered counter-clockwise) of HH 202, HH 529, HH 203, 
HH 204, HH 528, HH 269, and HH 625 intersect, this gen- 
eral area having been previous ly identified as the most 



likely source for these outflows (pally et al.|2000| |Q'Dell 



Doi 2003 Henney et al. 20071). This figure is similar 



JBalty 
flTTr" 



to Figure 5, except the color balance is slightly different, 
we have added labels identifying the individual sources, 
and tangential velocities are shown projected for only 100 
years. 




Fig. 7. — This 38.7 // x41.0 // sample of the central region shown 
in Figure 5 uses the same color coding (although a different color 
balance) and the identification numbers of the sources have been 
added. Motions projected over 100 years are now shown. The 
symbols have the same meaning as in Figure 3. 

Although we have derived higher accuracy tangential 
velocities, these have not narrowed down the identifica- 
tion of the sources of the large outflows except that vari- 
ous outflows seem to originate in the eastern part of this 
small field. If the features we call in this paper the HH 
269-Ram and HH 529-Ram are from the same source, its 
position is narrowed down to a Right Ascension interval 
of only 15". However, to link these two flows as a sin- 
gle bipolar flow demands that one of the sides has been 
deflected towards the observer by about 116°. Although 
a mechanism that provides deflection as a bea m passes 
through a d ensity gradient has been proposed ([Canto fc 



Raga 



19961), it seems a stretch to use this as the pro- 
cess that allows linking HH 269 and HH 529 into a single 
bipolar flow. The infrared source 146-351 and far infrared 
source 144-351 are the top candidates for the sources of 
these outflows. Although the other outflows must also 
arise in this same region, linking them all into a single 
source requires invoking the same large amount of deflec- 
tion multiple times, which is unlikely. Multiple infrared 
and radio sources exist in this region and if they indi- 
vidually produce bipolar flows, then the observational 
selection effect of seeing only blueshifted outflow from 
a source lying behind the main ionization front would 
apply and explain the dominance of blueshifts. If the ex- 
tinction that causes the sources to be visible only in the 
infrared was local to the source, then one would expect 
to see both redshifted and blueshifted ext ended flow. 

A rece nt infrared polarization study (|Hashimoto et 
ah] |2QQ7|) may provide additional useful! 



Hashimo to et al.| ( |2QQ7| ) argue that imbedded stars with 
bipolar outflows create dumbbell shaped cavities along 
the inner parts of the outflow and that starlight scattered 
from the surface of these cavities would have a charac- 
teristic polarization pattern. They have studied the in- 
frared polarization in this region and find polarization 
signatures around 136-355, 145-356, and 144-351, mak- 
ing them prime candidates for outflows. Possible outflow 
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Fig. 8. This 7.1" x 7.1" view of HH 208 shows the structure and 
changes in the object. The left two panels are brightness inverted 
views of the inner and outer regions of the object in the GO 10921 
F673N ([S II]) PC image (0.0455"/pixel). The right hand panel is 
a positive depiction of the ratio of earlier to recent images, showing 
the changes of structure in the inner region and general expansion 
in the outer. The nearby star is 137-222=JW 420( |Jones fcW alker 
1985) and the vertical axis points towards PA=45 . 

from 145-356 or the star 145-355 only 0.9" north of it 
is discussed in § 3.8. The direction of the putative flow 
from these two stars is indistinguishably the same as the 
high velocity flow feeding HH 203. The east side of the 
polarization pattern associated with 144-351 is generally 
oriented to be like a flow in the direction of the HH 529- 
Ram features, thus strengthening the argument that this 
is the source for the HH 529 flow. 

In the southwest portion of Figure 7 one finds the 
sources of the molecular outflows, 136-400 producing the 
bipolar CO flow, whose northwest blueshifted portion 
creates the optical object HH 625 plus the series of H2 
knots ahead of it, and 136-355, producing the SiO out- 
flow, portions of which may be in the features we are 
calling the western end of HH 269. The polarization pat- 
tern for 136-355 looks perpendicular to the SiO outflow 
direction, so its interpretation is not clear. Neither of 
the brighter radio sources (137-408, CS 3; 134-411, FIR 
4) in the corner of Figure 7 have clearly associated op- 
tical outflows, although the shocks constituting HH 530 
and shown in Figures 3-5, appear to move away from the 
region of 137-408 and 134-411. 

3.4. Outflow from the BN-KL Region 

Study of the BN-KL region is primarily the subject of 
infrared and radio investigations since most of outflow- 
ing material and all of the luminous central sources are 
invisible at optical wavelength. However, opti cal stud- 
ies have aided in narrowing down the origin (Jones & 
Walker|1988"l|Doi et al.|2"002| . Study of the motion of the 
three most massive sources in BN-KL sho ws that they 
are expandin g away from a c ommon point (Rodriguez et 
all|2005| IGomez et al.||2005|) which we label as "DYN- 



CTR" in our figures. The timescale fo r expansion of th e 



optically visible shocks is 1000 years (Doi et al. 2002), 



although thi s may be shorter if t he m aterial ha s been 



decelerated ( |Lee fc BurtonpOOOp . |G6mez et ah] ( |2005D 
find a dynamic timescale of about 500 years. The cause 
of this explosive decay and the physics relating it to the 
molecular outflow are still unresolved, but it is reason- 
able to link the outflow t o whatever occurred ther e. An 
alternative view is that of |Beuther fc Nissen| ( |2008| ), who 
argue that SMAl is the source of the high velocity flow, 
with a dynamic timescale of about 1000 years, and that 
source I is the source of the low velocity outflow. 

It is worth asking whether the optical features can help 
to illuminate what is happening. Many of the molecular 
fingers (H2 is the best way of tracing them) have low ion- 
ization optically visible shocks at their tips (HH 201, HH 
210, and the HH 205+HH 206+HH 207 combination are 




Fig. 9.— This 46.3"x 35. 1" WFPC2 F673N inverted brightness 
image at O.OQQG^/prxel from program GO 10921 shows the region 
near HH 208 with north up. DYN-CTR indicates the posi tion of 
the dynami cal center of expansion of the radio sources (Gomez 
|et al.| |2005) and the nearby sources B N, I, and n. The submm 
source bMAl's position is taken from Beut her &; Nissen| J 2008). 
The m otions of objects 126-214 and 129-216 are trom |Doi et al.| 
( 2002 ) and the others are from this study. The dashed line indicates 
the axis of symmetry of the shock 136-233. The light diagonal band 
in the upper left occurs at the juncture of two of the WFPC2's 
detectors. 

the best exa mples), with HH 201 being the br ightest and 
best studied ( [Graham et aLl|2003| [Doi et al ||2004[ ). The 
resolved form of the shocks reflect their oblique angle to 
the viewer, but one of them, HH 208, appears closest to 
the origin, is seen almost along the line of sight, and fell 
within the high resolution (0.0455"/ pixel)field of view 
in our GO 10921 observations, having previously been 
imaged in the lower resolution (0.0966"/ pixel ) WFPC2 



C CD's in programs GO 5469 and GO 8121 flO'Dell et 
al. 1997|) . It is undetectable against the nebula in the 
Ha, [JN II], and [O III] filters but stands in good con- 
trast in the [S II] and [O I] filters, especially [S II]. 
There is no discernible motion of this object, the cen- 
troid (5:35:13.59 —5:22:23.9) of this irregular object be- 
ing constant to about 3 km s _1 . What is most notewor- 
thy is the change of structure that has occurred. Figure 
8 shows the present appearance both in the core (middle 
panel) and the outer parts (left panel), while the right 
panel shows the ratio images between the recent (GO 
10921) and oldest (GO 5469) observations (12.76 years), 
after magnifying to the same resolution. The changes 
in structure are remarkable, well beyond what might be 
attributed to sampling with different size pixels, and ar- 
gue that the form of this nearly head-on bow shock is 
highly variable. The outer parts of it are expanding at 
about Vt = 48 km s _1 , yielding a dynamical timescale of 
70 years. We have determined the radial velocity o f the 
object from [S II] spectra ( [Garcia-Diaz et al~ 2008) and 
found Vq = —165 km s~\ high , but much lower than in 
HH 201 and HH 210 ( jHu|1996[ ). These numbers indicate 
that the motion is within 1° of our line of sight and it 
must be nearly in the same direction as its source. 

A wider field of view around the putative source is 
shown in Figure 9. There is a series of low ioniza- 
tion features moving away from the source region, two 
knots within them having measured tangential motions 
(timescale 900 years for motion from the DYN-CTR 
point). However, this region is very complex, there being 
a linear array of fine features passing from the proplyd 
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Fig. 10. — The central 7.9" x 7 '.!" region around the proplyd 
driving HH 502 is shown with PA=-23.4° as up. The left panel 
shows the object in inverted brightness and the right the ratio of 
the first-epoch/second-epoch images. The white lines indicate the 
measured motion away from the central star over the observation 
interval of 4.29 years. The upper left knots have a dynamic age of 
16±2 years and the lower right knots a dynamic age of 56±4 years. 
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Fig. 11. This inverted brightness, 60.9" x 47.6", image of 
LL6 and its driving jet (HH 876) in the ACS-WFC F658N filter 
shows the motions in the jet and the knots in the curved shock 
characteristic of the LL Ori class of objects. 



154-240 flSchultz et al.||1999] ), passing through the star 
image 147-234 and SIVlAl, south of DYN-CTR, and stop- 
ping at the bright knot of [S II] emission, short of reach- 
ing HH 208, but pointing towards the shocks containing 
129-216 and 126-214. There is a well defined shock (136- 
233) moving toward s PA=198°, and it ha s an H 2 feature 



on its leading edge (Schultz et al. 
axis does not point back to any of 



1999| ). Its symmetry 
he candidate sources. 



If it began at DYN-CTR, the timescale would be 700 
years at Vt=38 km s~ , so that again one finds a dy- 
namic timescale for the motion that is about the same as 
the outer HH objects and the central dynamical event. 
Repeated observations of this region with WFPC2 will 
not be possible since that instrument will be replaced 
during servicing mission SM4 (now sched uled for Octo- 
ber, 2 008), but repeat observations in H 2 (Schultz et al. 
1999) using the NICMOS instrument may be able to use 
tangential velocities to narrow down the source of the 
outflow. 

3.5. Outflows in the HH 502 Region 

The area south of the Huygens regi on is known to be 
very rich in bow shocks and outflows ( Bally fc Reipurth 
2Q01J [Bally et al.||2001| [2006] |HenrTey et al.|l2007| ). This 
is probably because the luminosity of a shock is pri- 
marily dependent upon the density and velocity of the 
outflowing material, which must be basically indepen- 
dent of distance from ^Ori C, whereas the confusing 



background brightness of the nebula has diminished be- 
cause of the larger distance from the source of ultravio- 
let photons that ultimately determine the emission-line 
brightness of the nebula. The regi on centere d on HH 502 
has been thoroughly studied by |Bally et aL] ( |2006[ ) using 
ACS-WFC F658N images taken 1.465 years apart. In 
the present investigation the images were made with a 
separation of 4.293 years, hence the tangential velocity 
resolution was superior, being 24.3 km s~ per 0.05" per 
pixel. The larger field of view over which we determined 
tangential velocities is shown in Figure 6 and smaller en- 
closed fields around the proplyd 280-931 and LL 6 (HH 
876) are shown in Fig ure 10 and 11. We adopt the sys- 
tem of nomenclature of Bally et al. (2006), which suffices 
as only a few new features were measured. In fact, the 
tangential velocities are typically high enough that our 
longer time-base observations produc e results (Appendix 



C) that are very similar to those in |Bally et al. (2006) 
with the exception of o bject HH 502 -S7, where our value 



is about half of that of Bally et al. (2006). 



HH 502. — The most striking system in the measured 
field is HH 502, the bipolar flow discovered fro m ground- 
based observations (Bally & Reipurth 2001), which in 
form very much resembles HH 513 in the Huygens region 
both in having quite broad shocks, but also in showing a 
slight curvature away from ^ 1 Ori C with increasing dis- 
tance from the originating proplyd 280-931. 

One observes the structured jets that drive both out- 
flows as one approaches 280-931. Figure 10 shows both 
the first epoch appearance of this inner region and also 
shows in a ratio image the changes in structure and move- 
ment that has occurred in even the relatively brief in- 
terval of the observations we have used. The dynamic 
timescale for these objects are 15.1 years (1), 17.2 years 
(2), 54.8 years (14), and 56.9 years (15), giving us the 
shortest time resolution of outflow events in the Orion 
Nebula. The timescales for the outer features are much 
longer, being N3 (4490 years), N2 (1790 years), Nl (87 
years), S3 (460 years), S4 (420 years), S5 (850 years), S6 
(1550 years), and S7 (3430 years). 

HH 876-LL 6. — This object closel y resembles the 



archetype for these objects LL Ori (Bally et al. 2000 
2006 Henney et al.| [2007) in having a broad bow shock 



that is indistinguishably stationary facing the flow of gas 
in its region, and a nearly perpendicular bipolar jet that 
produces its own set of rapidly moving small bow shocks. 
In this case the bow shocks have the short timescales of 
275 years (HH 876-3), 140 years (HH 876-2), and 56 years 
(HH 876-1). 

3.6. HH 540 

The "Beehive" proplyd (181-826) HH 540 and its sur- 
rounding field has been imaged with the ACS -WFC in 
the F 658N filter and thoroughly discussed in |Bally et al. 
(2005). There is an inner silhouette disk with orthogo- 
nal microjets associated with shocks along slightly curved 
trajectories to both the north and south. However, the 
most remarkable thing about the object is the series of 
rings along the bright shell of the proplyd, these being 
symmetric about the bipolar flow axis, suggesting the 
nickname Beehive as it resembles a beehive as rendered 
in old cartoons. The cause of these rings is unknown, 
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Fig. 12. This 8.8"xll.8" image of HH 540 made with the 
WFPC2 is color coded as described in the text. The vertical axis 
points toward PA=355.3°. The inner silhouette disk is discernible 
around the central star as are two microjets, the lower (southern 
and opposite to the ionizing star ^ 1 Ori C) being brighter. The 
ripples in the ionized shell are probably caused by irregularities in 
the rate of mass loss within the microjets. 



but Bally et al. (2005) have probably identified the most 
likely cause. The jet that penetrates the ionized portion 
of the proplyd must be irregular in flow, thus providing a 
variable stimulus for perturbation s that would the n pass 
along the local ionization front. |Bally et al.| ( |2005| ) point 
out t hat a single simila r ring is seen in the proplyd 182- 
413. Bally et al. (2005) calculate propagation velocities 
for the perturbations from the dynamic ages of the dis- 
tant shocks driven by 181-826's outflow, finding velocities 
of a few kilometers per second and linking this to sound 
speed in the gas just inside the ionization front of the pro- 
plyd. It seems equally likely that the rings one now sees 
are produced by the much closer (and younger) knotty 
structure in the outflow, which would yield higher veloc- 
ities of the rings and arguing that they are perturbations 
propagating in the ionized portion of the proplyd. 

Our WFPC2 images cannot improve on the clarity of 
the view of HH 540, however, our observations included 
images in a wide range of emission line filters (F502N, 
F656, F658N, and F673N) in addition to the medium- 
width filter F547M . This allows us to discuss the ion- 
ization structure better, especially since the ACS-WFC 
F658N used in the earlier study passes both Ha and [N II] 
emission equally well, thus leaving an ambiguity about 
what one is actually observing. We have created a color 
image (Figure 12) showing F656N (Ha) as blue, F658N 
[N II] as green, and F673N ([S II]) as red. There is the 
unusual situation that the micro jet is more visible on the 
side away from 1 Ori C. This is especially true in [S II] 
emission, which indicates that this side is shielded by the 
inner disk from the FUV radiation tha t can pass through 
the proplyd's ionization front, whereas pally et aL (2005) 
explain the greater brightness in F658N on the south side 
as being due to there being less obscuration by the disk 
on that side due to the orientation. It is more likely that 
the EUV illumination process is more important as the 
south side is the blueshifted. One also sees in the color 
image that there seems to be a progression of ionization 
states in the rings. 

We have used the superior emission-line isolation to 
calculate the surface brightness of the proplyd in Ha, 
usin g the technique and calibration constants presented 
in |Q Dell fcT)oi| ( |l999| >. The full chord width is 3.5", 
the distance from the central star to the bright cusp is 
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Fig. 13. This series of 5.6"x 9.3" WFPC2 images of HH 510 
shows the wide variation in appearance in different emission line 
filters. PA=46° is up. The arrowhead indicates motion over 20 
years. 

1.0", and the distance from the brigh t cusp to the end 
of the tail opposite ^Ori C is 6.3". |Bally et al | ( |2005 ) 
indicate that the proplyd is 303" fr om its photoiomzmg 
star ^Ori C. Using the definition of |Q' Dell fc Doi| j 1999 ) 
for the average surface brightness (their § 3, paragraph 
3), we find this to be 8.3 x 10 8 Ha photons cm 



2 s- 1 



ster . A dir ect comparison with the calibration in Bally 



et al.| (12005), where the Ha line could not be isolated 
indicates that their surface brightness is about 1.7 times 
too small. 

We can compare the results for this large distant pro- 
plyd with similar objects found closer to ^Ori C by refer- 
ence to a study of the Huygens region proplyds ( jO'Dell 
1998). The chord size fits right into the middle of the 



pattern of size versus distance. The surface brightness is 
0.14 times (dex=0.85) that predicted by the upper limit 
in the surface brightness-distance relation. This relation 
would only be satisfied if the proplyds all lie in the same 
plane as 1 Ori C and all of the ionizing Lyman contin- 
uum photons were absorbed at the ionization front of the 
proplyd. The former condition is not satisfied, causing 
the Huygens region extinction corrected proplyds to run 
about 0.35 dex lower than the maximum brightness pre- 
diction. Reconciling the observed surface brightness and 
that expected from the properties of the Huygens region 
proplyds requires an extinction of 0.50 dex at Ha. We 
do not know the extinction in this region, but a value 
of 0.50 dex from the foreground veil is similar to that in 
the outer regions wh ere it has been determined (O'Dell 
fc Yusef-Zadeh|2000 ). The extinction would be even less 
if the proplyd lies more than an average amount out- 
side the plane containing ^Ori C. These comparisons 
indicate that the Beehive proplyd quantitatively behaves 
like the inner region proplyds in spite of its much greater 
distance from ^ 1 Ori C and its bizarre appearance. 

3.7. New, Expanded, and Newly Clarified Herbig Haro 

Objects 

The availability of this new set of data on tangen- 
tial motions and the high resolution ACS-WFC images 
has allowed the discovery of several new outflows, under- 
standing new detail in a few HH objects, and has caused 
rethinking about what composes two already designated 
objects. 

HH510.— T his proplyd (109-327) was discovered ( |Bally| 



et al. 2000) to have a microjet on a wide bandpass 
(F606YV) VVFPC2 image and designated as HH 510. The 
short microjet extended in the direction of two sharp fea- 
tures (si and s2 in Bally et al.'s Figure 15 and Table 3). 
Our new observations were able to measure motions in 
the brighter s2 feature and to determine that it is at the 
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Fig. 14. 



This inverse brightness depict ion in the ACS- 



WFC F658N (Ha+ [N II], JHenney et aL][2007t ) filter shows an 
86.9" x 66.3" field near the head of HH 528(large dark feature in 
the middle) with north up. The insets show enlarged 10 /; x 10 /x 
regions. The proplyd 1 57-533 and its mi crojet has previously been 
designated as HH 512 ( |Bally et al.|2 000). The alignment with this 
microjet, their forms, and their motions all indicate that the shocks 
179-515 and 206-454 are part of this HH object. 

tip of the high ionization jet (the s2 feature is visible in 
Ha, [N II], and [O HI]). This knotty jet can be traced all 
of the way (5.9") to s2 (Figure 13). Feature s2 is moving 
at 101 km s~ towards PA=64°. The timescale for this 
motion is only 120 years and presumably much shorter 
for the knots closer to the proplyd. 

HH 512.— T he proplyd 157-53 3 (Figure 14) was previ- 
ously known (Bally et al. 2000) to have a microjet pro- 
jecting towards PA= 60°, which lead to its designation 
as HH 512. Examination of the GO 10246 ACS-WFC 
F658N image shows there to be a very small (0.9" length) 
inner silhouette disk indistinguishably perpendicular to 
the microjet and that the microjet extends to 3. 5" from 
the central star. A shock (17 9-515) within th e "cap" 
of HH 528 was already known (Doi et al. 2002) to pos- 
sess an anomalous motion (we now measure this to be 
104 km s~ towards PA=56°), but a connection with HH 
512 was not made. We now note that this bow shock lies 
along the projection of 157-533's microjet and is moving 
away from that proplyd. Our reprocessing of the images 
made near HH 203 and HH 204 show an additional bow 
shock (206-454) lying almost on the line between 157-533 
and 179-515 and is again moving away from the source 
(103 km s~ towards PA=76°) and is almost certainly a 
component of this HH object. The timescales are compa- 
rable to those of many other Huygens region HH objects, 
being 770 years (206-454) and 1700 years (179-515). The 
remarkable thing about these bow shocks is that their 
small size is evidence that the driving jet, which is es- 
sentially unresolved near the source, has not expanded 
to more than 1.4" at a distance of 84.3". 



H H 513 .— This well studied (Bally et al. 2000 Doi et 



aL|2002 ) (Figure 3, Figure 4, and Figure 5) object is the 
only bipolar object in the Huygens region where one can 
measure tangential motions on both sides of the outflow. 
Our observations indicate that the east side components 
are moving at about 34 km s~ and those on the west 
side at 75 km s~ , indicating timescales of 100, 200, 240, 
310, 430, 1500, and 2700 years for the seven component 
shocks. We have measured [O III] spectra for the west 




Fig. 15. This 49.7 // x23.2 // image of the components of HH 
514 and HH 518 uses the same WFPC2 images as Figure 5 but 
with a greyscale rendering and has PA=314.4° as up. The bright 
star just off the field of view is ^ 1 Ori C. The length of the arrows 
indicate tangential motions in 100 years. 



flow from the dataset created by Doi et al. (2004) and 
determined that Vq = —30 km s~ , indicating that the 
west side is the closer side and that it is at an angle of 
about 37° with respect to the plane of the sky. 



HH 514 and HH 518.— These join HH 873 ([Garcia-Diaz 



Henney 2007) as the only redshifted Herbig Haro out- 



fiows in the Huygens Region. The features included in 
these two HH objects have been differ ent in the several 
articles that have considered them ( Bally et al.| |2000 



sly, 



our 



Doi et al.||2002| |Q'Dell fc Doil|2003[ ). Fortunate 
new data allow clarification of this rather complex re 
gion, which is shown in Figure 15. The simplest system 
is HH 514, which begins as a redshifted microjet aris- 
ing from the proplyd 170-337 (there is a faint blueshifted 
counter jet in the opposite direction) and continues with 
a group of rapidly moving knots ( 170-334) about 3" north 
of the proplyd ( |Bally et al.|2000[ ). |Bally et al |fl2000[ ) then 
say that there are a series of associated shock features at 
5.6", 8.9", 12.9", and 14.6". Only the feature at 5.6" 
(which is visible on our images but its motion was not 
measured) is probably part of HH 514 unless ther e is a 
very abru pt change in direction of the flow. Ba lly et| 



al 



( 2000| ) define HH 518 as a single shock (designated 
Here as 163-342) associated by symmetry and proximity 
with the proplyd 159-350 and say that there were other 
fainter associated shocks furthe r from the sourc e, with- 
out making clear identifications. Doi et al. (2002) argued 
that the shock at 172-326 (172-327 here) was a new HH 
object, calling it HH 605. However, with only little cur- 
vature, it should be i dentified as part of the HH 518 flow 



(O'Dell & Doi 2003). We believe that HH 518 is com- 



posed of the lour shocks 177-322, 172-327, 165-334, and 
163-342 for which we have determined tangential veloc- 
ities. The higher ionization radial velocity mapping of 
this region (Doi et al. 2004[ ) provide radial velocities of 
V Q = 98 km s _1 , 104 km s _i , 74 km s~\ and 102 km s" 1 
respectively. This means that HH 518 is also anomalous 
within the Huygens region in being a redshifted f low. 



Moreover, e xamination of the Ha spectrum 159-350 (Doi 
et al. 2004) shows a redshifted jet. The values of de 



rived from our tangential velocities and these radial ve- 
locities are —51°, —69°, —55°, and —77°, that is, the flow 
is highly inclined and pointing away from the observer. 
A very similar result applies for HH 514, where our tan- 
gent iaj^noti£n^nd_the radial velocity (^0=152 km s~ ) 
from |Ballyet al.| ( |2000D give 9 = -69°. The unusual and 
similar values for these flows argues that there is some- 
thing different and possibly related in the outflows from 
159-350 and 170-337. 
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itions in 100 years are shown. 



Fig. 16. This 30.1" x 23.7" ACS-WFC F658N image (north 
up) shows the newly designated HH 995 and HH 997. HH 997 
originates in the dark-disk proplyd 143-522, where one sees a faint 
microjet extending to the northeast, perpendicular to the dark 
disk and towards the two associated shocks. The middle shock 
(149-513) is centered on a narrow linear feature that is likely to be 
an extension of the microjet. The microjet forming HH 995 arises 
from 141-520. 

HH 626. — This object is an incomplete double ellipse 
(possibly two opposing systems o f broad bow s hocks) 
centered on the proplyd 143-425 flO'Dell fc Doi] [2003). 
Our tangential velocities indicate timescales of 450 and 
500 years for the northern components and 360 and 750 
years for the southern components (excluding a single 
northern component of 1900 years, which is probably due 



to an error in the Vt= 8 km s determination). 

HH 995 and HH 997. — Our ratio images indicate a 
pair of shocks (152-508 and 149-513) oriented towards 
PA=45° and moving towards PA=51°. The westward of 
the pair has a three arcsecond linear feature lying along 
its symmetry axis (Figure 16). A projection backwards 
through these two shocks passes close to the newly dis- 
covered bright proplyd 143-522 with a silhouette dark 
disk, which shows a microjet extending (PA=45°) to- 
wards the shocks and is almost certainly their driving 
source. There is another newly discovered similar object 
(141-520) with a microjet extending towards PA=40° 
which has been designated as HH 995. The timescales 
for movement from 143-522 are 1100 years (152-508) and 
300 years (149-513). 

HH 998.— The bright rim proplyd 160-353 (Figure 17) 
reveals a new feature on the GO 10921 and the ACS- 
WFC images in that it has a 0.4" microjet extending to- 
wards PA=47°and nearly opposite to this microjet there 
is the rapidly moving (80 km s _1 ) shock 158-355. It is 
not unusual in Orion to find that the microjet is on the 
side facing the nearby photoionizing star (^Ori C in this 
case) and there is no reason to not expect a similar, but 
less visible, outflow in the opposite direction. It is this 
currently-invisible counter-jet that must be driving 158- 
355 and the timescale for the separation is only 100 years. 

3.8. Possible new HH Objects 

There are many features for which we have determined 
tangential velocities but these features have not been as- 
signed to existing or new (HH 995, HH 997, and HH 998) 
HH objects. In some cases there are indications of the 
need to assign them membership in a new HH object, 




Fig. 17. This 8.65" x 10.35" ACS-WFC F658N image (north 
up) shows the newly discovered microjet extending to the north- 
east from the proplyd 160-353 and the shock 158-355 which is in 
the opposite direction of the microjet. This ensemble is newly des- 
ignated as HH 998. HH 529-1 is the first of a series of large, high 
ionization, eastward moving shocks within the system designated 
as HH 529. The bright feature extending to the northeast from 
159-350 is an overexposed image of a companion star, not another 
microjet. The arrow designates the motion in fifty years. 

but the level of evidence does not reach a sufficient level. 
We mention several of these in this section. 

Outflow from 117-352? — Two nearly aligned shocks (110- 
352 and 116-352) are moving away from the nearby pro- 
plyd 117-352 (Figure 5), however, this area overlaps with 
the large HH 269 shocks and they may be part of that 
HH object. 

Outflow from 150-355?— There are two shocks (149.3- 
351.7 and 150.5-351.9) just north and moving away from 
the star 150-355 (Figure 7). This motion is quite different 
from other flows occurring in this area. 

Outflow associated with the Dark Arc? — There are four 
low tangential velocity features (130-345, 131-342, 133- 
342, and 140-344) which appear (Figure 7) to be moving 
away from the geometric center o f this unexplained ob- 
ject QO'Dell fc Yusef-ZadehpOOOl ). 



Motions near 143-353 and 145-355+145-356. — There are 
three southeast moving features (151-359, 154-401, and 
158-405) aligned with both the infrared stars 145-355 and 
145-356 and the more distant radio source 143-353). A 
northwest moving object (144-354) lies between the two 
possible driving sources. If all four objects are part of 
the same bipolar flow, then the source would be either 
145-355 or 145-356. In the infrared polarization study of 
Hashimoto et al. (2007) one sees a polarization signature 



around the latter pair of stars. Although they favored 
145-356 as the source, the low angular resolution cannot 
make a definitive distinction between the two. The di- 
rection of the southeast moving features points towards 
the high velocity jet that drives HH 203 and it is possible 
that these are all one system. 

A Very High Velocity Feature. — One of the highest tan- 
gential velocity features found in this study is 137-317 
(Vt = 155 km s~ , PA=113°). It is very low contrast, 
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only about 10% brighter than the surrounding back- 
ground and is only clearly discernible in F658N, barely 
detected in F631N and F673N, and not visible in F502N 
and F656N. It is a very low ionization shock of about 
1.5" width, symmetric to the direction of motion. There 
are no candidate sources nearby for this unusual feature. 

3.9. General Properties of Stellar Outflows 

The particular advantage of studying stellar outflows 
from objects in the ONC is that one is dealing with a 
group of stars all about the same age and sharing similar 
histories. This means that we can draw general conclu- 
sions that probably are applicable to other HH, which 
are usually found as isolated, single, objects. 

A striking pattern of the HH objects in the ONC is 
the intermittent nature of the outflows, a not uncommon 
chara cteristic that w as first documented on s maller scales 
by | Reipurth | ( |1989[ ) and on larger scales by Reipurth et 
L997J). Our time resolution is probably best for the 
shocks and jets because distant bow shocks can ac- 
tually be the result of multiple close-together in time 
outflows. In single objects like HH 502 we find well de- 
fined variations in the bipolar outflow with time gaps of a 
few tens of years. These staccato intervals of outflow are 
probably spaced out at larger intervals, producing the 
larger scale features that one sees at greater distances 
from the sources. This means that bipolar outflow is 
not at all continuous, rather, there is evidence that it 
occurs in concentrated periods at intervals of a few hun- 
dred years and during the period of outflow there is a 
irregular mass loss with a timescale of about ten years. 
Even shorter timescale variations may be present, but 
our observations would not have been able to detect typ- 
ical outflows of shorter than a few years. 

With the exception of HH 625, the largest scale HH ob- 
jects in the Huygens region appear to originate in a small 




area, immediately northeast of the brightest Orion- S ra- 
dio and infrared sources and they do not arise fr o m the 
SiO and CO molecular outflows. Zapata et al. (2006) 
point out that the sources driving the known molecular 
outflows are not particularly luminous and it is not sur- 
prising that even less luminous and less obscured sources 
would be able to produce the highly visible objects like 
HH 202, HH 203, HH 204, HH 269, and HH 529. The 
ubiquity of blueshifts in these HH objects argues that 
these originate from different young stars, because link- 
ing them to a common source would require large changes 
in orientation in the plane of the sky and along the line- 
of-sight with timescales of a few hundred years. 

These observations do not give us definitive informa- 
tion about long timescale variation in the direction of the 
axis of outflows since the visible outflow material is mov- 
ing through the ambient gas of the nebula (which will 
have a local low velocity flow as it moves away from the 
main ionization front and from the brightest part of the 
nebula) and, in some cases, will be affected by the stellar 
wind from the hot massive stars. 
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APPENDIX 
A. THE DISTANCE TO THE ORION NEBULA 

The exact distance to the Orion Nebula Cluster eludes determination because i ts distance is too great to simply 
apply traditional parallax methods and it is ill-advised to apply indirect methods (Muench et al. 2008 O'Dell et al. 
2008). The pre- main Sequence stars are not suitable because of the many variables in determining their properties 
and the main sequence O and B stars are the closest exemplars, so that neither group can be used for an accurate 
spectroscopic parallax. Jefferies (2007) has determined a distance using the statistical properties of the rotation of the 
cluster stars, but this result is very model dependent. 

Trigonometric parallaxes are th e fundamen tal method for determining distances, but even the Hipparchos satellite 
optical results are not definitive . Bertout et al.| (|1999|) noted one member star (HD 37061) and found a distance of 



1+168 



Wilson et al. 



(b 



(2005) used multiple stars in this direction and by grouping according to if they are 



361I^ 8 pc, while 

reddened by OMC-1 found 465+'^ pc. 
A rec ent optical interfer ometric study of 6 )1 Ori C derived a distance as a by-product of the study of this star as a 



binary (Kraus et al. 



They favored the Or 



2007). They found two equally good solutions, Orbit 1 (434±12 pc) and Orbit 2 (387±11 pc). 



3it 1 result, but only because of better agreement with previous determinations. 



Radio interferometric parallaxes are more promising. |Genzel et al.|p9 81) applied what is in-effect the moving cluster 
parallax metho d to H2O masers an d derived a value of 480±80 pc. More recently, radio trigonometric parallaxes have 



been reported. Hirota et al. (2007) have studied a H2O mas er in BN-KL, dete rmining its position in two dimensions 
from 15 observations and deriving a distance of 437±19 pc. |Sandstrom et al. (2007) measured a non-thermal source 
(GMR-A) in the inner ONC and found a distance of 389^ 2 i P c - However, their source is highly variable in flux, is 
partially resolved, (it also shows variations in its structure), and both of the reference sources are resolved. In a ddition, 
the radial v elocity of GMR-A is Vq = 14 ±5 km s _1 , which is quite differ ent from the value (V q = 25.6 km s~ Fiiresz 
et al.|[2008|) for the stars in the cluster. The most recent study is that of Menten et al. (2007) who measured only the 



Right Ascension component of the parallax of four non-thermal sources (including GMR-A). The observations were 
made at only four epochs and one of the sources was not detected during one epoch, so that parallaxes were determined 



HIGH VELOCITY FEATURES IN THE ORION NEBULA 



17 



TABLE 4 
Tangential Velocities in the Huygens Region 







V T 


PA 








V T 


PA 




Designation 


HH 


(km s- 1 ) 


(°) 


Filters 


Designation 


HH 


(km s - 


X ) (°) 


Filters 


100-246 


202 


49 


302 


F502N,F656N 


170-359 


529-III 


41 


119 


F631N,F656N,F658N,F673N 


104-245 


202 


38 


298 


F502N,F656N 


109-321 


625 


25 


302 


F631N,F656N,F658N,F673N 


112-246 


202 


31 


310 


F502N,F656N 


110-323 


625 


39 


307 


F631N,F673N 


115-255 


202 


57 


308 


F502N 


111-324 


625 


36 


289 


F631N 


116-250 


202 


58 


336 


F502N,F656N 


112-324 


625 


37 


285 


F656N,F658N,F673N 


116-252 


202 


36 


343 


F502N 


140-419 


626 


6 


001 


F656N 


116-255 


202 


80 


338 


F656N,F658N,F673N 


141-419 


626 


29 


001 


F502N,F656N 


116-256 


202 


41 


305 


F502N 


141-421 


626 


19 


343 


F502N,F656N 


116-254 


202 


52 


334 


F656N,F658N 


142-428 


626 


10 


163 


F656N 


117-249 


202 


36 


322 


F658N,F673N 


143-429 


626 


26 


208 


F502N 


117-255 


202 


33 


316 


F631N 


172-320 c 


726 


98 


094 


F502N 


116-344 


269 


52 


265 


F502N 


149-513 


997 


70 


045 


F656N,F658N 


116-346 


269 


38 


268 


F502N,F656N 


152-508 


997 


37 


057 


F656N,F658N 


118-345 


269 


61 


288 


F502N,F656N 


158-355 


998 


80 


242 


F502N 


125-346 


269 


38 


291 


F656N 


110-352 d 


NEW(?) 


33 


302 


F658N,F673N 


137-351 a 


269 


16 


283 


F502N,F656N 


116-352 d 


NEW(?) 


42 


269 


F631N,F658N,F673N 


138-350 a 


269 


33 


272 


F656N,F673N 


149.3-351. 7 e 


NEW(?) 


26 


334 


F656N 


139-351 a 


269 


11 


302 


F631N,F656N, F658N,F673N 


150.5-351.9 e 


NEW(?) 


27 


354 


F502N,F656N 


113-325 b 


510 


101 


64 


F656,F502N,F658N 


130-345 


Dark Arc 


20 


358 


F502N,F656N 


179-515 


512 


104 


056 


F656N,F658N 


131-342 


Dark Arc 


13 


316 


F502N,F656N 


206-454 


512 


103 


076 


F656N,F658N 


133-342 


Dark Arc 


20 


295 


F502N 


159-237 


513 


59 


277 


F502N,F658N 


140-344 


Dark Arc 


26 


026 


F502N 


161-236 


513 


60 


283 


F658N 


144-354 f 


NEW(?) 


50 


307 


F673N 


162-316 


513 


98 


273 


F502N,F658N 


151-359 f 


NEW(?) 


19 


109 


F673N 


168-235 


513 


47 


076 


F502N,F656N,F658N 


154-401 f 


NEW(?) 


26 


154 


F673N 


171-234 


513 


49 


053 


F502N,F658N 


158-405 f 


NEW(?) 


35 


116 


F631N,F656N,F658N,F673N 


173-233 


513 


18 


035 


F502N,F658N 


115-359 




32 


278 


F656N,F658N 


177-229 


513 


15 


025 


F502N,F658N 


129-237 




85 


351 


F502N 


170-334 


514 


48 


357 


F502,F656N,F658N 


131-246 




9 


226 


F658 


137-508 


516 


33 


142 


F631N,F673N 


132-343 




17 


46 


F502N 


163-342 


518 


17 


046 


F656N 


133-344 




46 


1 


F656N 


165-334 


518 


34 


035 


F502N,F656N 


134-332 




23 


268 


F631N,F658N,F673N 


171-328 


518 


33 


060 


F502N 


134-344 




34 


1 


F502N 


172-327 


518 


27 


056 


F502N,F656N,F658N 


135-243 




29 


318 


F502N 


177-322 


518 


58 


080 


F502N 


135-347 




22 


5 


F502N,F658N 


182-337 


523 


74 


096 


F502N,F631N,F656N,F658N 


136-233S 




38 


198 


F673N 


165-456 


528 


20 


190 


F631N,F673N 


137-317 




155 


113 


F658N 


166-434 


528 


19 


196 


F631N,F673N 


115-359 




32 


278 


F656N,F658N 


166-440 


528 


26 


164 


F631N,F673N 


143-326 




36 


289 


F502N 


166-442 


528 


28 


203 


F656N,F658N 


153-305 




77 


334 


F502N 


169-446 


528 


20 


171 


F656N,F658N 


155-253 




13 


13 


F502N,F658N 


170-459 


528 


19 
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F673N 


155-257 




39 


334 


F502N,F658N 


171-502 


528 


20 


161 


F631N,F656N,F658N,F673N 


157-326 h 




41 


305 


F502N 


177-454 


528 


37 


199 


F673N 


161-325 




11 


247 


F502N,F658N 


181-510 


528 


26 


154 


F631N,F673N 


161-410 




26 


208 


F658N 


184-513 


528 


18 


163 


F631N,F673N 


165-317 




25 


316 


F502N 


149.1-352.9 


529 


85 


104 


F502N,F656N,F658N,F673N 


165-406 




48 


96 


F656N,F658N 


149.6-353.4 


529 


95 


105 


F673N 


166-308 




98 


42 


F502N 


150.1-353.3 


529 


91 


105 


F631N,F656N,F658N,F673 


167-300 




23 


38 


F502N 


151.1-353.9 


529 


63 


90 


F502N,F631N,F656N 


169-325 




30 


125 


F502N 


152-354 


529 


18 


64 


F673N 


171-218 




17 


5 


F502N,F658 


155-354 


529 


18 


91 


F502N 


178-246 




19 


359 


F658N 


156-356 


529 


45 


109 


F502N,F656N 


178-246 




15 


6 


F658N 


161-353 


529-1 


30 


80 


F656N 


178-306 




20 


200 


F502N 


161-354 


529-1 


78 


100 


F502N,F631N,F658N,F673N 


180-245 




18 


23 


F658N 


167-358 


529-11 


92 


91 


F502N 












Note. — Obj 


ects are 


grouped by assigned Herbig Haro number or potential nev 


j grouping. When there is no 


candidate 


Herbig H 


'aro flow designation, they are listed 


individually. Wit 


hin each 
b Featui 


cateeorv obiects are ordered in increasing Right Ascension 


jet? Driven by proplyd 117- 


352? e D 


riven by ] 




Irregular Jet? 


•e S2 of|Ball> 


' et al.| 


|2000| c Outflow from LV2's (167-317) 


L50-355? 


f Driven by 144 


1-355? 8 


■Driven by BN(?), c 


.f. text. h Symmetric with LV6 (158-326) 











for only three sources, finding 417±7 pc. If GMR-A is excluded, then the remaining two sources give 412±6 pc. 

We have derived a distance of 436±20 pc by assigning weight 3 to the Hirota et al. 2007 study, weight 2 to the 
Menten et al. (2007) and Kraus et al. (2007) studies (Orbit 1), and weight 1 to the Genzel et al. (1981) studies. We 
adopt the rounded value of 440 pc in this paper. 



B. TANGENTIAL VELOCITIES IN THE HUYGENS REGION 



Insert Table [4] here. 
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TABLE 5 
Tangential Velocities near HH 502 



Bally Designation Tangential Velocity (km s 1 ) Position Angle (°) 



Comment 



HH 876-1 (west) 


117 


HH 876-1 (middle) 


67 


HH 876-2 


64 


HH 876-3 


53 


HH 560S-1 


14 


HH 560S-2 


33 


HH 560S-3 


48 


HH 502-1 


248 


HH 502-2 


329 


HH 502-N2(8) 


56 


HH 502-N2(7) 


79 


HH 502-N2(6) 


106 


HH 502-N3(10-ll-12) 


30 


HH 502-N3(NE) 


48 


HH 502-Nl(5) 


76 


HH 502-Nl(3) 


224 


HH 502-Nl(4) 


182 


HH 502-S3 


107 


HH 502-14 


102 


HH 502-15 


153 


HH 502-S4(16) 


158 


HH 502-S5(17) 


122 


HH 502-S6(19) 


100 


HH 502-S7(21) 


54 


HH 874(22) 


123 


HH 874(23) 


154 


HH 874(24) 


168 


HH 874(25) 


55 


HH 874(26) 


58 



274 
270 
260 
205 
303 
310 
300 
27 
27 
63 
118 
98 
24 
22 
43 
35 
35 
198 
207 
207 
202 
199 
197 
196 
126 
116 
157 
165 
164 



Divided into two components. 
Divided into two components. 



Additional feature measured. 



Note. — The Bally Designations are those given in Table 3 and Table 4 of |Bally et al.| {2006). 
C. TANGENTIAL VELOCITIES NEAR HH 502 
Insert Tabled here. 
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